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Abstract 

This paper presents a construction of isospectral problems on the torus. The con- 
struction starts from an SU(n) version of the XYZ Gaudin model recently studied 
by Kuroki and Takebe in the context of a twisted WZW model. In the classical 
limit, the quantum Hamiltonians of the generalized Gaudin model turn into classical 
Hamiltonians with a natural r-matrix structure. These Hamiltonians are used to 
build a non-autonomous multi-time Hamiltonian system, which is eventually shown 
to be an isomonodromic problem on the torus. This isomonodromic problem can 
also be reproduced from an elliptic analogue of the KZ equation for the twisted 
WZW model. Finally, a geometric interpretation of this isomonodromic problem is 
discussed in the language of a moduli space of meromorphic connections. 
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1 Introduction 



It has been argued for the last several years that isomonodromic problems are closely 
related to the Knizhnik-Zamolodchikov (KZ) equation One of the earliest observations 
is due to Reshetikhin |2| . Reshetikhin considered the Schlesinger equation, and concluded 
that the KZ equation may be viewed as a quantization of the Schlesinger equation. Harnad 
reformulated Reshetikhin's observation in a Heisenberg picture |§ . Let us recall here that 
the Schlesinger equation is an isomonodromic problem on the Riemann sphere. Since the 
KZ equation on the Riemann sphere can be generalized to the Knizhnik-Zamolodchikov- 
Bernard (KZB) equation M on the torus, one will naturally expect that an associated 
isomonodromic problem should exist on the torus. Korotkin and Samtleben, indeed, 
derived such an isomonodromic problem from the KZB equation of the SU(2) WZW 
model |J. Recently, Lavin and Olshanetsky proposed a general framework for this type 
of isomonodromic problems on a general compact Riemann surface . 

We present below an SU(n) and "twisted" version of the isomonodromic problem of 
Korotkin and Samtleben. The word "twisted" means that our isomonodromic problem 
is related to the "twisted WZW model" recently studied by Kuroki and Takebe [0]. The 
method of construction, too, is considerably different. Korotkin and Samtleben start from 
a Hamiltonian formulation of the Chern-Simons theory (related to the ordinary untwisted 
WZW model), and formulate the isomonodromic problem in terms of a meromorphic con- 
nection induced on the torus. We follow a more direct approach that has been known for 
a class of isomonodromic problems on the Riemann sphere |8], [|. This class of isomon- 
odromic problems can be systematically derived from a class of isospectral problems as 
a non-autonomous analogue; of particular interest is the case where the isospectral prob- 
lems are the so called Hitchin systems |10| and their generalizations to punctured Rie- 



mann surfaces [lT|. As such an isospectral problem, we now take an SU(n) version of 
the XYZ Gaudin model considered by Kuroki and Takebe in the study of the twisted 



WZW model. Like the ordinary Gaudin (or "Calogero-Gaudin" ) models fl2| , |T3 , H] , this 
generalized Gaudin model is a quantized Hitchin system in the sense of Beilinson and 



Drinfeld ||15|| . In the classical limit, mutually commutative Gaudin Hamiltonians turn 
into Poisson-commutative Hamiltonians with an r-matrix structure. We construct from 
these Hamiltonians a non-autonomous multi-time Hamiltonian system, then rewrite it 
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into a Lax system of isomonodromic type, and finally confirm that this isomonodromic 
problem is linked, in the sense of Reshetikhin, with Etingof 's elliptic analogue of the KZ 
equation (T^j (the "elliptic KZ equation" in the terminology of Kuroki and Takebe [[?]]). 
As shown by Kuroki and Takebe, this equation plays the role of the KZ equation for the 
twisted WZW model. 

Isomonodromic problems on the torus (and more general compact Riemann surfaces) 
have also been studied by Okamoto |L7|], Iwasaki ]nj and Kawai |19| by complex analytic 
and geometric methods. Although their work has been mostly focussed on scalar Fuchsian 
equations, their methods can be applied to matrix systems like ours. We shall show that 
this reveals a geometric origin of the Hamiltonian structure of our isomonodromic problem. 



2 Generalized Gaudin Model 

We now briefly review Kuroki and Takebe's generalization of the XYZ Gaudin model to 
an SU(n) spin system 0. 

Let us first introduce basic functions and matrices. Let X be the torus (elliptic curve) 
with modulus r, i.e., X = C/Z + Zr. For integer indices (a, b), let 9[ a b](z) denote the 
special theta functions 

0[ab]{z) = 9a _i i_b(z,r), 

7i 2 ' 2 n 

9 KK '(z, t) = exp[niT(m + k) 2 + 2Tci(m + k)(z + «')]. (1) 

Furthermore, let J a b and J ab be the n x n matrices 

J ab = g a h\ J ah = -J~ b \ (2) 

n 

where 

g = diag(l,uj,uj 2 ,---,u n - 1 ), h=(S i ^ j ), uj = e 2 ^ l ' n . (3) 

Note that g and h obey the algebraic relation gh = uhg. It is well known that these J's 
give a basis of su(n) (over R) and sl(n, C) (over C). J ab and J ab are "dual" bases in the 
sense that Ti(J ab J cd ) = 5 ac 6 bd . 

These functions and matrices are building blocks of Belavin's Z n -symmetric i?-matrix 
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and the associated r-matrix [20]. The .R-matrix is given by 



R(\)= £ W ab (X,i])J ab <g> J ab , (4) 

(a,6)6Z„xZ„ 

where 

W«{X,,) = (5) 
e [ab]{V) 

These Boltzman weights are n-periodic in a and b, thereby the summation is over Z„ x Z n . 
The r-matrix is the leading nontrivial part in the //-expansion of -R(A) at r) — 0: 

r(A)= Yl w ab {\)J ab ®J a \ (6) 

(0,6)^(0,0) 

where 

^ ) [afe](0)6 l [oo](A) ' 

and the prime means A-derivative, ' = d/dX. The summation in (a, b) is now over Z n x 
Z n \ {(0,0)}. We shall frequently omit showing this range explicitly (or, equivalently, 
obey the convention that woo(A) = 0). The r-matrix satisfies the classical Yang-Baxter 
equation 

[r ( 13 ) (A) , (p)] = - [ r ( 12 ) (A - fj), r (13) (A) + r ( 23 ) (//)]. (8) 

Here, as usual, the superscript means in which part of the tensor product C n <g> C n ® C n 
the r-matrix acts nontrivially, e.g., r^ 12 ^(A) = yj w ab (\) J ab (g> J a6 ® J. 

The generalized Gaudin model is a limit, as rj — > 0, of an inhomogeneous SU(n) spin 
chain with iV lattice sites. The monodromy matrix of this inhomogeneous spin chain is 
given by 

T(X) = L N (X-t N )---L 1 (X-t 1 ), (9) 

where t^s are inhomogeneity parameters (which eventually play the role of time variables 
in our isomonodromic problem), and Lj's are L-operators of the form 

L t (X) = £ W ab (X) J ab ® Pt (J ab ) (10) 

(ab) 

that act on the tensor product C n (g> Vi of C n and the representation space of an 
irreducible representation (pi,Vi) of su(n). These L-operators satisfy the well known 
equation of ll RLL = LLR" type with R being the above .R-matrix. The monodromy 
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matrix thereby acts on C n <g> V (V — V\ <S> • • • <S> Vjv). The leading nontrivial part in the 
77-expansion of T(A) at i] = can be written 

N 

T W = E E w *b(\ - u)j ab ® Pt (j ab ), (11) 

i=l (ab) 

which satisfies the fundamental commutation relation 

[T« (A) , T< 2 > (//)] = -[r(A - /i), T« ( A) + T (2) (/,)]. (12) 

The superscript now stands for the component of the two C ra 's in C n <g)C™<g)V in which the 
monodromy matrix acts nontrivially, e.g., T^(A) = X)i Z)( a &) ~ ii)^ab <8> / <E> Pi(J ab )- 

Mutually commutative Hamiltonians Hi (i = 0, 1, • • • , N) of the generalized Gaudin 
model are defined by the relation 

1 N N 

- Tr T(A) 2 = ]T CiV(X - U) + ^C(A - U) + H , (13) 

2 nxn — * — * 

t=l t=l 

where V(z) and ((z) are the Weierstrass functions with modulus r. Note that Cj's are 
quadratic Casimir elements of the algebra generated by pi(J ah ): 

Ci = \Y.Pi( J «>)Pi( jah )- (14) 
The ifj's for i = 1, • • • , N can be written rather explicitly, 

H i = E E™<*& - tjMJ^PjiJ*), (15) 

but explicit expressions of become far complicated. There is a linear relation among 
these Hamiltonians: X^Li #i = 0. 

3 Construction of Isomonodromic Problem on Torus 
3.1 Classical limit and Poisson bracket 

We now introduce the following matrix as a classical analogue of the T-operator of the 
generalized Gaudin model: 

N 

M W = EE W ^( A _ ti)J ab Af. (16) 

i=l (ab) 
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Af b, s are scalar functions of ti, • • • , ijv and r. Our isomonodromic problem will be formu- 
lated as differential equations for these functions. (Note that Af"s depend on t whereas 
pi{J ab Ys do not. As we shall show later, this amounts to the difference of the Heisen- 
berg and Schrodinger pictures in quantum mechanics.) We apply the same rule of raising 
and lowering the indices as the J-matrices to these coefficients. Thus, e.g., A,^ ab are 
defined by A^ ah = nui ab A^ a ~ b in accordance with the transformation rule for the J's, 

T „, ,ab j—a,—b 

Ji,ab — nUJ J 

The commutation relations of pi(J ab ) , s can be now translated into Poisson commuta- 
tion relations among the A-coefficients. This Poisson structure can be packed into the 
well known relation 

{M(A) f M(n)} = -[r(A - /i), M(A) ® J + J <g> M(fj)]. (17) 

The left hand side, as usual, is an abbreviation of E{M ab {\), M cd (fi)}J ab <g> J cd , where 
M ab are the coefficients of the expansion M = J2 M ab J^. In terms of the residue matrix 

A t = Res M(\)=J2J ab A?, (18) 

X ~ U (ab) 

the Poisson structure is nothing but the one induced from the Kirillov-Kostant Poisson 
structure on sl(n, C). Symplectic leaves of this Poisson structure are the direct product 
Oi x • • • x On of coadjoint orbits 0, in s\(n, C) on which Ai is living. These symplectic 
leaves become the phase spaces of the following non-autonomous Hamiltonian system. 

3.2 Non- autonomous Hamiltonian system 

Hamiltonians of the classical Gaudin model can be obtained as follows: 

i N N 

- Tr M(A) 2 = ]T dV(X - ti) + £ HiC(X - U) + H . (19) 

i=l i=l 

Cj's are Casimir elements in the above Poisson algebra. The Hamiltonians for i = 
1, • • • , N, as in the quantum case, can be written explicitly: 

H i = E E w ab(U ~ tj)A hab Af. (20) 

j&i) (06) 

Hq is also a quadratic form of Af's, though its explicit expression is complicated. 
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It should be noted that t^s in these formulas are just parameters, not time variables, 
of the classical Gaudin system. The classical Gaudin system is a Hitchin system on 
the punctured torus X \ {ti, ■ ■ ■ , tjv} |13|, |14| , thus ti's are nothing but the position of 



punctures. We now consider a non-autonomous multi-time Hamiltonian system with the 
same Hamiltonians, which gives a non-autonomous analogue of this Hitchin system. 

The non-autonomous Hamiltonian system possesses the puncture coordinates ti, ■ ■ ■ , t N 
and the modulus r as time variables. The equations with respect to the puncture coordi- 
nates are given by 

dA ab 

-^- = {Af,H i }. (21) 

The equation with respect to the modulus, in contrast, turns out to take the somewhat 
strange form 

dA ab N 

= {A?, (H - m 5>J7i)/W=T}, (22) 

where rjx is the constant (depending only on r) that arises in the transformation law 

C(z+l) = C(z)+m (23) 

of the Weierstrass ( function. 
3.3 Lax representation 

In order to confirm that the above non-autonomous Hamiltonian system is an isomon- 
odromic problem, we now rewrite it into a Lax form. 
A key is the following relation: 

1 N 

{M(A),-TrMGu) 2 } = ^^-AW/J-^J^MW]. (24) 

Z i=l (06) 

This relation can be reduced to a collection of cubic relations among w a bS, and can be 
proven by the same function-theoretic method as the proof of the classical Yang-Baxter 
equation (which reduces to quadratic relations among w a bs). 

The above formula can be used to rewrite the right hand side of the Hamiltonian 
equations into a matrix commutator. First, the residue of both hand sides of the above 
formula at // = ti immediately gives the relation 

{M(\),H t } = -[M\),M(\)}, (25) 



where 

A(X) = Y,w ab (X-t i )J ab Af. (26) 

(ab) 

(We have used the reflection property W- a - b (— A) = —w ab (X) of wWs, too.) The 
Poisson bracket with H requires lengthy calculations, which eventually boil down to the 
relation 

N 

{M(\),H }=\4TrV^lB(\)-r ll Y,t l A l (\),M(\)}, (27) 



i=l 



where 



N 

B{X)=Y.H Z ab^-U)JabA i t h \ 

i=l (ab) 

7 m _ w ab {\) ( 0' [ab] (X) gjaq(0) \ , . 

Za6(A) " W^T ivKA) " «o) ) ■ (28) 

The above formula for {M(A),Mo} is interesting in itself from two points of view. 
Firstly, this formula can be rewritten 

N 

{m(A),^ -^E^} = [47rv /z T5(A),M(A)], (29) 
i=i 

thus the previous strange linear combination of H^s in the Hamiltonian equation with 
respect to r emerges quite naturally. Secondly, the functions Z ab are also basic constituents 
of the elliptic KZ equation of Etingof (see Section 4). 

By this formula, the previous non-autonomous Hamiltonian system can be rewritten 
into the following Lax equations: 

^) = _ [A( A),M(A)]-^), 

»jjw =&BW , MW]+2 °m. ( 30) 

The second terms on the right hand side originate in differentiating w ab (X — tj) in M(A) 
by ti and r. The r-derivatives are converted into A-derivatives by the heat equation 

de [ab] (X) 1 d 2 6 [ab] (\) 

dr 47rv^T 8\ 2 ' 1 j 
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3.4 Isomonodromic Property 

As usual, the above Lax equations are integrability conditions of a linear system: 

(^-M(A))y(A)=0, (J- + ^(A))r(A)=0, (|- - 2S(A))y(A) = 0. (32) 

The first equation is a matrix ODE on the torus with regular singular points at ti, ■ ■ ■ , tjv- 
The other equations may be interpreted as isomonodromic deformations of this ODE in 
the following sense. 

Firstly, the solution of the ODE transforms along a small closed path around X — ti 

as: 

Y(X) - K(A)IV (33) 

Ti represents the local monodromy around A = £j. Similarly, the solution of the ODE 
transforms along the fundamental cycles a : z — > z + 1 and (3 : z —>■ z + r, respectively, as: 

F(A + 1) = hr x Y{\)Y a , 

Y(\ + r) = gY(X)T p . (34) 

The extra prefactors h -1 and g originate in the monodromy of M(A) along those cycles, 
i.e., M(A + 1) = h- 1 M(X)h and M(A + r) = gM(X)g~ 1 . T a and r> represent the global 
monodromy. "Isomonodromic deformations" means that these monodromy data are left 
invariant under deformations. 



4 Relation to Elliptic KZ Equation of Twisted WZ W 
Model 



The elliptic KZ equation of Etingof [T(J can be written 
d_ 



d 

(«7aT+ E EM*<-*i)M^)Pi(^ a6 ))-P'(*ir--,*iv) = 0, 

i(^) (06) 



a N 

E ^Z o6 (t i -i> i (J a6 )p i (J o6 ))F(t 1 ,---,i JV )=0, (35) 

" T i,i=l (a6) 

where k = + n, k is the level of the twisted WZ W model, and the above equation 
characterizes iV-point conformal blocks with the irreducible representations px,---,p N 
sitting at the marked points ti, ■ ■ ■ , t^. 

9 



Following Reshetikhin 0, We now add one more marked point at A with the funda- 
mental representation (C n ,id). The associated N + 1-point elliptic KZ equation becomes 

r) N 

(n— + £E^&( A - U)JabPi(J ab ))G(\,t 1 , ■■■,t N ) = 0, 

i=l (ab) 
j&i) (ab) 

+ J2 W ab(t t ~ X) Pl (Jab)J ab )G(X,t 1 , --•,%) = 0, 

(ab) 

d_ v 

dr 



i,j=l (aft) 

iV AT 

+ E E M*< - A)p t ( Jaft) J a " + E E - WabPAJ *) 

i=l (aft) i=l (aft) 

+ X; ^aft(O) J a ft J a6 )G(A, ttf) = 0. (36) 

(aft) 

Let us now take an invertible "fundamental solution" F = exp 5 § of the Appoint 
elliptic KZ equation, and consider equations for X = F~ 1 G. This is a kind of "gauge 
transformation" to the AT + 1-point elliptic KZ equation, and the outcome is the equations 

a N 

(«^ + E E W ab(\ ~ U)J ab A?)X(\, tjv) = 0, 

i=l (aft) 

J- - E«W A " *i) JobAf )X(X,t l , ■■■,t N ) = 0, 

1 (ab) 
Q N 

(«7T + 2 E E ^aft(A - ti) JabAf + X ^aft(O) J a ft J ah )X(A, * a , • • • , t N ) = 0, (37) 

° T i=l (aft) (aft) 

where 

Af = F(t x , ■ • ■ , t^pii^Ffa, •■■,%). (38) 

We thus obtain almost the same equation as the isomonodromic linear system in the last 
section, but notice that there are a few discrepancies: 

1. The last equation contains extra terms (a linear combination of J a bJ ab )- 

2. The Af b, s are not scalar functions, but operators on V = V\ ® • • • ® Vn- 

3. The above linear system contains an arbirary parameter k. The previous isomon- 
odromic linear system can be reproduced by putting k = — 1. 
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The first discrepancy can be readily remedied, because the extra terms are scalar 
(J a bJ ab — I l n ) an d can be "gauded away" by a scalar gauge transformation X i— > e^^X 
with a function /(r) of r only. This does not affect the other part of the above equations. 

The second discrepancy is more essential. From the point of view of Reshetikhin 
and other people || [f| , the above system is a quantization of the isomonodromic problem 
in the last section. (This is parallel to the relation between quantum and classical Hitchin 
systems.) The operators A® b depend on t's and r, but inherit the same commutation 
relations as the J-matrices from pi(J ab ). The passage from pi(J ab ) to Af 3 amounts to 
the change from the Schrodinger picture to the Heisenberg picture. "Classical limit" 
now means replacing these operators by functions on a phase space. This gives our 
isomonodromic problem. 

The parameter k plays a role in the limit towards the critical level k = —n. In 
this limit, the (classical or quantum) isomonodromic problem in the above sense turns 
into an isospectral problem, i.e., a Hitchin system on the punctured Riemann surface 
X \ {ii, • • • , t/v}- A more careful analysis of this limit leads to the so called "Whitham 
dynamics" of the spectral curve, which describes the isomonodromic problem as a slowly 
varying isospectral problem || 21 . 



5 Geometric Origin of Hamiltonian Structure 

Our isomonodromic problem may be reformulated in a geometric language at least in two 
different ways. One option is Levin and Olshanetsky's framework |J based on Hamilto- 



nian reduction. Another option is Iwasaki's geometric framework ||18|| ) which has been 
successfully applied by Kawai to a scalar isomonodromic problem on the torus. In 
this section, we follow the second approach and show a geometric origin of the somewhat 
involved Hamiltonian structure of our isomonodromic problem. 

In the geometric approach of Iwasaki and Kawai, meromorphic linear ODE's on a 
compact Riemann surface X are converted into meromorphic connections V on a holo- 
morphic vector bundle E. All singular points are assumed to be regular singular points. 
One can then develop an analogue of the Kodaira-Spencer theory on the moduli space M. 
of meromorphic connections with a given constant set of local monodromy data around 
the regular singular points. The tangent spaces of this moduli space, i.e., infinitesimal 
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deformations leaving invariant the local monodromy, can be described by the language 
of twisted de Rham cohomologies with coefficients in the endomorphism bundle End-E. 
A key of Iwasaki's ideas is that the Poincare-Lefschetz duality in those twisted de Rham 
cohomologies induces a closed 2-form (or a symplectic form) Q on the moduli space M.. 
This closed 2-form describes the Hamiltonian structure of isomonodromic deformations. 
Kawai generalized these ideas of Iwasaki to a setting where the Riemann surface itself 
also deforms, and considered, as an example, the case of isomonodromic deformations of 
a second order Fuchsian ODE d 2 y/d\ 2 = Q(\)y on the torus. It is this generalization of 
Iwasaki's framework by Kawai that we now attempt to apply to our problem. 

It is rather straightforward to reformulate our problem into Kawai 's framework. In 
our case, the endomorphism bundle EndE in twisted de Rham cohomologies is replaced 
by the Lie algebra bundle sl(n, C) tr of Kuroki and Takebe 0. This Lie algebra bundle is 
obtained from the trivial bundle sl(n, C) x X over the universal covering X = C of X by 
identifying (A, A) ~ (hAhr 1 , A + 1) ~ (g _1 Ag, A + r). The closed 2-form Q is defined by 
the Poincare-Lefschetz pairing < 5\M, 5 2 M > of two infinitesimal deformations 5±M (A) 
and &2M(A) of M(A) that leave invariant the local monodromy around £1, - • • , t^: 

£l(5 1 M, 5 2 M) =< 8 1 M, 5 2 M > . (39) 

Calculation of this pairing, as Kawai illustrated, can be reduced to a kind of residue 
calculus, and eventually yields the following expression of Q: 

N N j N 

n = Y / dH l A dtt + d(H - Vl J2 UHi) A - = - J2 Tr dB i A dA i- ( 4 o) 
i=i i=i 27TV-1 i= i 

Here A4 is the residue matrix of M(A) at A = ti, and dBi is linked with dAi as 

[Ai, dBi] = dA,^ (41) 

This result is almost parallel to Kawai 's result, but the third part on the right hand 
side does not exist in the case of Kawai's isomonodromic problem for second order scalar 
ODE's. This part is nothing but the Kirillov-Kostant symplectic form on the symplectic 
leaves (coadjoint orbits) of sl(n, C), on which the residue matrix Ai moves as a function of 
time variables. A similar expression of the fundamental 2-form for general isomonodromic 
problems is derived in Levin and Olshanetsky's work |J in the framework of Hamiltonian 
reduction. 
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The above expression of f2 also clearly shows the correspondence between the Hamil- 
tonians and the time variables. In particular, one can see why we have to take the strange 
linear combination (H —J2iLi UHi) /2ny/— T of ifj's in order to describe the isomonodromic 
deformations with respect to the modulus r. 



6 Conclusion 

The point of departure of our construction is the generalized Gaudin model of Kuroki and 
Takebe. Its classical limit provides a Poisson commutative set of Hamiltonians with an 
r-matrix structure. The isomonodromic problem is first formulated as a non-autonomous 
multi-time Hamiltonian system, then rewritten into a Lax representation. From the Lax 
representation, we have been able to see a link with the elliptic KZ equation of Etingof, 
which advocates Reshetikhin's observation in a non-zero genus case. We have also found 
a geometric interpretation of the Hamiltonian structure of this isomonodromic problem 
in the framework of Iwasaki and Kawai. 

Our isomonodromic problem is a generalization of the Schlesinger equation on the 
Riemann sphere. It should be stressed that constructing such an isomonodromic problem 
on the torus is by no means an easy task. This seems to be one of reasons that complex 
analytic and geometric studies on isomonodromic problems have been mostly focussed 



on scalar equations fll7|, [18| , |19|| . The approach from the KZ equation is rather suited to 
dealing with matrix equations, and deserves to be pursued further. 

An interesting issue in this respect is to construct difference and g-difference analogues 



of isomonodromic problems from difference |22| and g-difference p3|, [24[] analogues of 
the KZ equation. For instance, Jimbo and Sakai constructed a g-difference analogue of 
the sixth Painleve equation by a g-analogue of the standard method for isomonodromic 
problems |[25|| . It seems likely that their g-difference equation is linked with a g-difference 
KZ (qKZ) equation. 

The author thanks Hironobu Kimura and Takashi Takebe for useful comments and 
discussion. 
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